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ABSTRACT 

We analyse a high-resolution, fully cosmological, hydrodynamical disc galaxy simula- 
tion, to study the source of the double-exponential light profiles seen in many stellar 
discs, and the effects of stellar radial migration upon the spatio-temporal evolution of 
both the disc age and metallicity distributions. We find a "break" in the pure expo- 
nential stellar surface brightness profile, and trace its origin to a sharp decrease in the 
star formation per unit surface area, itself produced by a decrease in the gas volume 
density due to a warping of the gas disc. Star formation in the disc continues well 
beyond the break. We find that the break is more pronounced in bluer wavebands. By 
contrast, we find little or no break in the mass density profile. This is, in part, due to 
the net radial migration of stars towards the external parts of the disc. Beyond the 
break radius, we find that ^60% of the resident stars migrated from the inner disc, 
while ^25% formed in situ. Our simulated galaxy also has a minimum in the age pro- 
file at the break radius but, in disagreement with some previous studies, migration is 
not the main mechanism producing this shape. In our simulation, the disc metallicity 
gradient flattens with time, consistent with an "inside-out" formation scenario. We do 
not find any difference in the intensity or the position of the break with inclination, 
suggesting that perhaps the differences found in empirical studies are driven by dust 
extinction. 

Key words: galaxies: spiral- galaxies: stellar content - galaxies: kinematic and 
dynamic - galaxies: structure - galaxies: abundances methods: numerical 



1 INTRODUCTION 



Since 1940, it has been known that the stellar disc 
light in spiral galaxies decrease s exponentially with radius 
(|Pattersonl 1 1940l : iFreemanll 19701 ) . This is thought to be the 
result of the initial angular momentum distribution of the 
gas c l oud that collapsed to form the disc llFall fc Efstathioul 
198d: iDalcanton et al.lll997l; iMo et al.|["l99Sl; Ivan den~Bosch 



20011 : iFerguson fc Clarkell200ll h However, since the pioneer- 
ing work of van der Kruit (1979), it has become clear that 
this exponential profile does not extend t o arbitrarily large 
radii in the majority of disc galaxies. iPohlen fc Truiillol 
(2006) and Erwin et al. (2008) showed - for later and ear- 
lier types, respectively - that, in fact, the galaxies for which 
this happens represent a minority, identifying three classes 
of surface brightness behaviour in the outskirts of disc galax- 
ies: ~60% of the galaxies show a deficit of light with respect 
to an extrapolated singular exponential (Type II); ~30% 
show an excess of light in the outskirts respect to a singu- 
lar exponential l|Erwin et al.ll2005l : IPohlen fc Truiillol |2006| . 
Ty p e III); while only ~10% show pure exponential pro- 
files l|Bland-Hawthorn et all 120051 : IPohlen fc Truiillol 120061 : 



lErwin et al1l2008l . Type I) out to >10 scale-lengthfS Fur- 
ther departure from the notion of the "disc-as-exponential" 
paradigm leads to even more elaborate classifications based 
upon, for example, the position of the "break" with re- 
spect to the outer Lindblad Resonance of the bar (for 
Type II) or the exact shape of the s urface brightness profile 
in the outskirts (f or Type III) fee.. IPohlen fc Truiillofl2006l ; 
lErwin et al]|200g| h 

The break radius (-Rbr), or the radius at which the 
light profile starts to deviate from a pure exponential is 
empirically seen to occu r at ~1.5— 4.5 disc scale-lengths 
l|Pohlen fc Truiillol l2006h . Outside the break, the surface 
brightness does not immediatel y drop to zero, but follows a 
secon d (different) exponential (Ide Griis et al.ll200ll ; IPohlen! 
2002). The fact that such breaks are very common, at least 
in late-type galaxies (Beckman et al. 2006), indicates that 
they either form easily or survive for significant periods of 
time. The exact causes though, 30 years on from the iden- 
tification of this "break" phenomenon, remain unclear. The 



1 These percentages are very dependent upon the morphological 
class of the galaxy (see Beckman et al. 2006). 
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leading candidates can be divided broadly into those related 
to angular momentum conservation and those related to a 
star formation threshold. 

In the case of the former, it can be shown that if angu- 
lar momentum redistribution does not occur in the disc, the 
outermost stellar radius will reflect the maximum value of 
the angular m omentum of the bar yonic material (that of the 
protogalaxy) ([van der Kruidll987f ). However, numerical sim- 
ulations have shown that non-axisymmetric instabilities do 
drive a substantial redistributi on of angular momentum in 
the disc (jPebattista et af . 2006) and that, in fact, this can be 
a responsible agent for a break. This is because the angular 
momentum redistribution leads to an increase in both the 
central densities and disc scale-lengths. As the angular mo- 
mentum redistribution cannot act efficiently to arbitrarily 
large radius, a break in the surface mass density distribu- 
tion results. 

Alternatively, breaks may be due to the effect of star 
formation thresholds (Kennicutt 1989). Within this latter 
scenario, the break radius would be located where the den- 
sity of the gas is lower than the critical valu e for star for- 
mation. The existence of ext ended UV discs (jThilker et al.l 
120071 ; iGil de Paz et al.l I2OO7T ), the lack of co rrelation be- 
tween Ha "cutoffs" and optical disc breaks (jPohlen et al.l 
120041 ; lElmegreen fe Hunter! [2006). and the exponential de- 
cay of the light outside the break, however, have all 
been used as arguments a gainst this picture. That said, 
lElmegreen fc Hunter! l|2006l ) showed that a double exponen- 
tial profile may result from a multicomponent star forma- 
tion prescription, where turbulent compression in the outer 
parts of the disc can allow for cloud formation and star for- 
mation despite sub-critical densities. More recently, models 
have shown that redistribution of stellar material in the disc 
can also lead to exp onential profiles out side the star forming 
disc's break radius ijRoskar et al.ll20o"8l ). 

Not surprisingly, hybrid models have been suggested 
which combin e aspe cts of both pictures. For example, 
Ivan den Boschl i|200ll ) investigated a combination of the col- 
lapse and the threshold model, claiming that gas breaks are 
a direct consequence of the angular momentum of the pro- 
togalaxy while the stellar breaks are determined by a star 
formation threshold. 

While each proposed model appears capable of ex- 
plaining the presence of disc breaks, it is equally im- 
portant that any successful hypothesis also takes into 

account the fact that: (i) breaks a re not seen in all 

galaxies ifrBarton fc Thompson! 1 19971 ; IWeiner et all l200ll ; 
Blan d-Hawthorn et al.l 20051 ); (ii) breaks are seen in discs of 
all morphological types, from SO l|Erwin et al.ll200"8h to Sm 
(eg. , Kregel et al. 2002; PT06) and even in irregular galax- 
ies l|Hunter fc Elmegreenl 12009) - albeit, not with the same 
frequency; and (iii) breaks are observed out to redshift z~ 1 
(|Perej|2004l : iTruiillo fc Pohlenll2005l ; lAzzollini et aill2008bl ) . 
Recent observational studies have now imposed additional 
constraints on the models for break formation. For exam- 
ple, de Jong et al. (2007) have shown that the break posi- 
tion is independent of stellar age, a result compatible with 
that of Bakos et al. (2008), who found the position of the 
break to be independent of the photometric band. It has also 
been found that the prominence or intensity of the break 
decreases with increasing distance from the galactic mid- 



plane l|de Jong et al. 2007 ) and for older stellar populations 
l|Bakos et al.ll2008l ; lAzzollini et al.ll2008al). 

Furthermore, the recent work of llloskar et al l (|2008l ) 
has demonstrated the importance of stellar migrations in 
defining the final mass density and surface-brightness den- 
sity profiles of disc galaxies. The re-organization of stars in 
the disc can have important consequences for the observa- 
tional constrains used for the chemical evolution models, as 
the age-metallicity relation or the metallic ity distribution , 
as previously shown by several studies (e.g. iHaywoodl 120061 . 
120081 ; ISchoenrich fc Binnevlliooll iRoskar et al.ll2008l) . 

Previous work in this area has concentrated upon some- 
what idealised scenarios of isolated galaxies. The advantage 
of this approach is that higher resolution can be achieved 
more readily, with a greater associated exploration of param- 
eter space. However, such an approach necessarily carries 
with it the disadvantage of not being able to take into ac- 
count various processes such as the cosmological infall of gas 
and/or the interactions with other galaxies that are intrinsic 
to the hierarchical assembly paradigm. Indeed, interactions 
with satellites can have an important influe nce in shaping 
the disc, as suggested by sim ulations (e.g., lYounger et alj 
120071 ; iPenarrubia et alj [2006) and should not be ignored. 
Our work aims to fill this hole by addressing the possi- 
ble physical origins of the two-component surface brightness 
profiles, but now within a cosmological context. We will ex- 
amine the likelihood that such profiles are the consequence 
of dynamical processes, star formation processes, or a com- 
bination of the two. 



2 SIMULATION DETAILS 

Our simulated disc was realised using multi-resolved, cos- 
mological simulations generated with the N-body, hydrody- 
namical code RAMSES ( Tevssierj f2002) which is based upon 
an Adaptive Mesh Refinement (AMR) technique. RAMSES 
includes gravitation, hydrodynamics, radiative cooling, and 
heating processes, immersed within a uniform, Haardt & 
Madau (1996) ionising radiation field. In addition to solv- 
ing the Euler equations with a net cooling term, RAMSES 
includes a phenomenological treatment of star formation, 
which contributes directly to the chemical enrichment of the 
interstellar and intergalactic medium via subsequent Type II 
supernovae0 The spatio-temporal evolution of the gas metal 
abundance, as well as the redshift-dependent photoionising 
background, are consistently accounted for in the computa- 
tion of the net cooling rates of the gas (Courty & Teyssier 
2009, in preparation). This combined contrib ution has been 
fit u sing the CLOUDY photoionisation code (|Ferland et al.l 
I1998D . 

The star forma tion prescription w i thin RAMSES is de- 
scribed in detail by iDubois fc Teys sier (2008); only a brief 
summary is provided in what follows. Star formation is per- 
mitted in cells whose density is higher than a given thresh- 



2 The current implementation of chemistry within RAMSES is re- 
stricted to the global metal content (Z), under the assumption 
of the instantaneous recycling approximation; an extension to 
RAMSES which relaxes this approximation and includes the effects 
of Type la supernovae and asymptotic giant branch stars is cur- 
rently under development (Few et al. 2010, in preparation). 
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old, according to the following rate: p\ — —p/t*, where the 
star formation timescale is proportional to the local free-fall 
time, t* = to(p/ po) 1 ^ 2 ■ As in Dubois & Teyssier's work, we 
set this timescale to to — 8 Gyr, with an associated density 
threshold for star formation of po = 0.1 H cm" 3 . Kinetic 
feedback of supernovae energy to the surrounding interstel- 
lar medium is modeled with a blast-wave solution that es- 
sentially mimics the expansion of superbubbles whose radius 
is fixed to two cells. Thermal feedback is accounted for by 
using a polytropic equation of state in the high-density re- 
gions. The amount of gas then turned into stellar material 
in each eligible gas cell is m,(l + r]sN +r]w), where the mass 
locked in long-lived stars is m*, and m, (77s at + r/w) is the 
mass carried away by the blast-wave, with the wind driving 
parameter r\w — 1, and the mass fraction of stars recycled 
into supernovae ejecta i]sn = 0.1. Photometric properties of 
the resulting stellar populations were derived by assuming 
that each star particle was a single stellar population (SSP), 
and then applying the models of Bruzual & Chariot (2003) 
with a Salpeter (1955) initial mass function (IMF) and the 
Bertelli et al. (1994) isochrones. 

We run our simulation to z = in a ACDM Uni- 
verse with the following cosmological parameters, fiA=0.7, 
Q. m = 0.3, n b =0.045, (78 =0.9 and H =73 kms _1 s _1 Mpc _1 . 
To carry out cosmological simulations of disc galaxies at the 
sub-kpc scale, a multi-resolved approach was adopted. Can- 
didate dark matter halos were chosen in a low-resolution, 
large-scale structure simulation. The size of this ACDM sim- 
ulation was 20/i _1 Mpc, contains 128 3 dark matter parti- 
cles, and was generated with RAMSES as part of the Horizon 
Projecl|f]. Our multi-resolved simulation used three nested 
boxes whose initial conditions were centred on our candi- 
date halo and whose resolution ranged from 128 3 in the 
external part of the computational volume to 512 3 for the 
most nested box. The latter resolution is also the one of 
the coarse grid in the central area and the adaptive mesh 
nature of RAMSES was employed to refine a further seven lev- 
els, thereby reaching a spatial resolution of 435 pc in the 
central area by z=0. Dark matter particles have a mass of 
6xl0 6 Mq and the initial gas mass per cell was ~10 6 Mg. 

As noted in Gibson et al. (2008), the selection of the 
candidate halo was made without any preconceptions re- 
garding spin A or triaxiality T. Care was taken though to 
ensure that the final halo was not contaminated by low- 
resolution dark matter particles up to at least five virial 
radii. The disc galaxy discussed here sits within a low-spin 
(A=0.02), mildly oblate triaxial (T=0.32) halo of dynamical 
mass 7.6xlO n M©. For redshifts z<2.3, the rotational axis 
of the disc could be aligned readily using the angular mo- 
mentum vector of the gas cells; the analysis which follows is 
therefore restricted to z<2.3. Further technical details of the 
simulation and alignment process are forthcoming (Courty 
et al. 2009, in preparation). 



3 CHARACTERISTICS OF THE DISC 

Hydrodynamical simulations of galaxies within the canoni- 
cal hierarchical structure assembly framework are known to 

3 http://www.projet-horizon.fr 



have a number of difficulties in reproducing disc-dominated 
galaxies (e.g., Springel & Hernquist 2003; Abadi et al. 2003; 
Bailin et al. 2005; Okamoto et al 2005; Governato et al. 2007; 
Scannapieco et al. 2008a). A major problem is that baryons 
condense early and then transfer a significant fraction of 
their angular mo mentum to the dark m atter as the final 
galaxy assembles l|Navarro fc Whitell994l ). As a result, these 
galaxies contain a significant fraction of their final baryonic 
mass in a spheroidal-like component supported primarily by 
velocity anisotropy, with consequent bulge-to-disc ratios in 
excess of those encountered in late-type spirals today. 

To derive the bulge-to-disc ratio of our simulated 
galaxy, we categorise the st ars a s "disc" or "sphero i d" fo l- 
lowing I Abadi et all (|2003bl ) and IScannapieco et all l|2008h . 
We derive the distribution of the circularity parameter e 
(where e = jz/jtir, and j z is the z-component of the spe- 
cific angular momentum of each star, where the z-axis is 
the symmetry axis and and j c ir is the angular momentum 
expected for a circular orbit at the same radius, r - i.e., 
jcir = rVdt (?*))• Figure [2] shows the e distribution of all 
stellar particles within a sphere of radius of 20 kpc at z=0 
(black). We have marked, with different colours, those par- 
ticles inside a sphere of radius 1.5 kpc (red), and those with 
radius 3< r <15 kpc and within 3 kpc of the disc mid-plane 
(green). Considering a stellar disc composed of stars with 
r<15 kpc, |zj<3 kpc, and 0.8<e<1.2, and including only 
those stars rotating in the plane of the disc (with cos a > 
0.7, where a is the angle between the angular momentum 
vector of the particle and the z-axis, after Scannappieco et 
al. 2008b), the disc-to-total ratio in our galaxy - measured 
as the ratio between the stellar mass on the disc and the to- 
tal stellar mass - is 0.37. Put another way, the bulge-to-disc 
ratio is, not surprisingly (and in keeping with the aforemen- 
tioned particle-based disc simulations), large with respect 
to that observed in late-type spirals (Sb-Sd), being more 
consistent with the ratio expected for an early-type (Sa-SO) 
galaxy. Figure [1] shows a bulge-to-disc decomposition of our 
simulation, using the K-band light profile and employing an 
exponential law for the disc and a Sersic law for the bulge. 
The fit was done in a iterative manner, following MacArthur 
et al. (2003), resulting in a photometric bulge-to-disc ratio of 
log(B/D)=— 0.5, consistent with that expected for an early- 
type (Sa-SO) galaxy (e.g. morphological type T~3: Simien 
& de Vaucouleurs 1986). 

The stellar bulge associated with our simulated disc 
is partially supported by rotation (with V ro t/cr~0.5, where 
Vrot is the maximum rotational velocity). In the simulation 
here, much as for the Milky Way bulge (admittedly, perhaps 
just coincidentally) , V ro t=70 km s _1 . 

Figure [3] shows a density map of the stellar disc from 
two different perspectives (face-on and edge-on). The disc 
extends to ~20 kpc, with an exponential scale-length (mea- 
sured in the V-band) of ~3.2 kpc, consistent with values 
found for late-type spirals such as the Milky Way (e.g. Juric 
et al. 2008 suggest a scale-length for the Milky Way's disc 
of ~2.6 kpc). 

While the stellar disc extends to ~20 kpc, the bulk of 
our analysis is restricted to stars within a radius of 15 kpc. 
This is mainly due to the "law-of-diminishing-returns" in the 
outskirts where the number density of stellar particles de- 
creases dramatically and the root-mean-square (RMS) dis- 
persion of relevant azimuthally-averaged quantities begins 
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Figure 1. K-band surface brightness profile of our simulated disc, 
plotted out to ~3 disc scalelengths. The fitted profiles of the disc 
(exponential law, shown in purple) and the bulge (Sersic law, 
shown in dark red) are overlaid, along with the sum of the two 
(shown in orange). 
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Figure 6. Azimuthally-averaged surface brightness profiles in the 
SDSS g-band, integrating particles within 1, 3, and 5 kpc from 
the mid-plane, as indicated in the inset. Error bars represent the 
RMS dispersion between the mean values found in eight arbitrary 
octants of the disc. 
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E = jz/jcir 



Figure 2. Distribution of e = jz/jclt f° r an stars within a sphere 
of r< 20 kpc (black line), within a sphere of 1.5 kpc (red line), 
and all those with 3<r<15 kpc and distance from the mid-plane 
<|3j kpc (green line). 



to increase rapidly. Furthermore, considering only particles 
inside this 15 kpc radius means we also reduce the possible 
contamination from co-spatial stellar halo particles. 

We end by noting that all "error bars" plotted through- 
out the paper are derived from the RMS associated with the 
mean values from eight arbitrary octants of the disc. 



4 CHARACTERISTICS OF THE DISC BREAK 

The first two rows of Figure [4] show the evolution of 
the azimuthally-averaged surface brightness profile in the 
V— band and the total stellar surface density profile, re- 
spectively, at five different epochs from 2=2 to z~0. It can 
be seen that the disc has already developed a break by 
z ~ 1 in the surface brightness profile and maintained its 
presence through to z=0. We fit simple functional forms 
(l-i(R) = /io + 1.086 x R/h) to the inner and outer expo- 
nentials and consider that the break occurs at the position 
where the functions intersect. At redshift z = 0, the break 
radius is i?br~9.5 kpc (i.e., ~2.9 disc scale-lengths). The 
break intensity - measured as the angle between the inner 
and outer exponential profiles - changes with redshift, but 



not in any obvious monotonic/systematic manner. It can 
also be seen that the break is much shallower in the total 
stellar mass surface density distribution. In fact, in most 
snapshots it is almost non-existent, a point to which we re- 
turn in Section [7. II 

Figure[S]shows spatially-resolved 2-dimensional star for- 
mation rate (SFR) maps of our simulation at z~0. The 
break radius is denoted by the vertical arrows (edge-on 
view) and circle (face-on view); an immediate point of in- 
terest is that there remains measurable star formation out- 
side the break radius. This star formation is not axisym- 
metric, with star formation extending beyond 7?br on one 
side, in particular. Such star formation in the outer disc 
(be yond the break radius) has been observed in many galax- 



ies dFergus on 1998; I Ferguson et al . 1998l; lThilker" et al 1 120051 . 



12007 ; Gil de Paz et al. 2007 ). Before assessing the physical 
origins driving the formation of the break in the surface 
brightness profile, we first review the main observational 
characteristics, to determine if our light profile is a fair rep- 
resentation of those observed in real galaxies. 



4.1 Dependence of the break upon distance from 
the mid-plane 

In the case of NGC 4244, Ide Jong et"aH (|2007l ) show that 
the break occurs at the same radius independent of height 
above the mid-plane. They further demonstrate that the in- 
tensity of the break decreases with height above the plane. 
Figure [6] shows the surface brightness profiles derived from 
our simulation, in the SDSS g-band, for samples of stellar 
particles at different dist ances from the plane . The empiri- 
cal characteristics seen bv lde Jong et all (|2007l ) are similarly 
seen in the simulation - specifically, the break radius posi- 
tion does not change with distance from the mid-plane, but 
it gets shallower. 



4.2 Dependence of the profile upon photometric 
band 

Different photometric bands are sensitive to different stellar 
populations. While the main contributors to the luminosity 
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Figure 3. Stellar density projections of the disc (left: face-on; right: edge-on). A 40x40x kpc region is shown in both panels. The circle 
in the left panel and the arrows in the right panel indicate the position of the break radius (see Sec. |4j. 
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Figure 4. Evolution with redshift - as indicated in the inset to each panel — of the V-band surface brightness profile (first row), total 
stellar surface density profile (second row), star formation rate surface density profile (third row), and gas surface density profile (fourth 
row), all on logarithmic scales. The profiles were derived using all stars with r<15 kpc and |z| <3 kpc. Dashed lines represent fits to the 
inner and outer exponentials in the V-band surface brightness profile. The break radii correspond to the intersection point of these two 
exponentials, and are noted by the vertical dashed lines in each panel. 



in the blue bands are the young starsQ this contribution 
decrea ses for the redder bands. Interestingly, iBakos et alj 
(2008) have found recently that the position of the break 
is independent of the photometric bandpass employed to 



4 Other hot stellar populations, including blue horizontal branch 
stars and blue stragglers, may provide an additional contribution 
in some systems (see, e.g., Rose et al. 1985; Maraston & Thomas 
2000; Trager et al. 2005) 



obtain the surface brightness profile, although break is shal- 
lower in the redder bandpasses. These results suggest that 
the break is due mainly to a truncation in the profile of the 
young stars and that the profile of the old stars should not 
show the same departure from a pure exponential. 

Figure[7]shows the surface brightness profiles of the stel- 
lar disc in the same two photometric bands (g and r) used in 
the empirical study of Bakos et al. (2008) , in addition to the 
redder K-band. In our simulation, the break appears at the 
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Figure 5. Edge-on and face-on 2-dimensional projections of the present-day integrated (over the last 100 Myr) SFR for our cosmological 
disc. For the face-on projection, only the stars at a distance of 5 kpc from the plane are included, to avoid contamination from the 
companion. The total line-of-sight in the edge-on projection is 40 kpc. The vertical arrows (edge-on) and circle (face-on) correspond to 
the position of the break radius in the V-band surface brightness distribution. 



same position (~9.5 kpc) but, in agreement with the obser- 
vations, is also shallower in the redder band (Ks) than it is 
in the bluer one (g). Furthermore, Bakos et al. (2008) trans- 
form the surface brightness profiles into total stellar mass 
density profiles, using mass-to-light ratios infer r ed fro m the 
colours, following the prescriptions of lBell et al. | (III). This 
process is admittedly not free of uncertainties, but, at face 
value, their results suggest that the mass density profile does 
not show a break, in agreement with that found for our sim- 
ulated disc (see Fig. [4| . 

Figure[8]shows the stellar mass surface density distribu- 
tion in our simulated disc for stars of different ages. While 
the position of the break is the same for those stellar pop- 
ulations showing a break in the stellar density distribution, 
not all stars show this truncated density profile. In partic- 
ular, stars older than ~8 Gyr show, if anything, a slightly 
"upbending" profile. Thi s latter result ag rees with a number 
of observational studies: iDavideel (2003) studied the stellar 
populations in the outskirts of NGC 2403 and M33, show- 
ing that while the number of young main-sequence stars 
present a "downbend" with radius, the number density of 
the more evolved stars in the red- and asymptotic- giant 
br anch phases do not ; Similar results were also obtained 
bv lGalleti et~al1 (I2004TI for M33 Conversely, for the case of 
NGC 4244. Ide Jong et all (20071 ) claim a sharp break even 
when only red giant stars (old stars for a constant star for- 
mation history) are considered. This issue clearly needs re- 
solving with a larger sample of galaxies. 



4.3 Colour profiles 



Recently, Bakos et al.l (2008) , for the local universe, and 
Azzollini et al .1 ll2008bl) . out to redshift z~l.l (see also 
Jansen et alT ~i|2000)) have found that the colour gradi- 
ents of galaxies with Type-II truncations are different from 
those with Type-I or pure exponential profiles (recalling the 
nomenclature introduced in Section 1). In particular, galax- 
ies showing Type-II truncations have a "U-shaped" colour 
profile, with the minimum (bluest color) at the position of 
the break. 

Figure [9] shows our predicted colour profile at four dif- 



ferent epochs from z=l to 2=0. In agreement with the ob- 
servations, our simulated disc also shows a U-shaped profile 
with the bluest colours near the position of the break. 



4.4 Age and metallicity profiles 

Colours are difficult to interpret in terms of simple stellar 
populat ions due to th e well-known "age-metallicity degen- 
eracy (Worthevlll9"94l ); in principle, the colour profile could 
be an effect of age, metallicity, or a combination of both. Fig- 
ure [TD] shows the mass- weighted stellar age and metallicity 
gradient for the final output of our simulation. It can be seen 
the age profile presents a minimum at roughly the position 
of the break. However, the metallicity profile is remarkably 
smooth all the way out ~15 kpc, without any visible change 
at the break radius. Therefore, the U-shape of the colour 
profiles seen in previous section is due to an age effect, and 
not to a metallicity effect. 



4.5 Evolution of the break with redshift 

Several studies have now analysed the e mpirica l evo- 
lution of the break radius with redshift (iPerezI |2004| ; 
iTruiillo fc Pohlenll2005l ; lAzzollini et al.l2008bl ) . Using a sam- 
ple in excess of 200 galaxies, for example, Azzollini et al. 
suggest that, for a given stellar mass, the radial position 
of the break has increased with cosmic time by a factor of 
1.3T0.1 between z ~ 1 and ~0. They also found that, in the 
same period of time, the evolution of the surface brightness 
level in the rest-frame B-band at the break radius (/ibr) has 
decreased by 3.3T0.2 mag arcsec -2 . 

Because the mass of our simulated galaxy evolves with 
redshift, we cannot compare directly the evolution of these 
parameters for a galaxy with a given mass. For this reason 
we re-analysed these data (kindly provided by the authors). 
We separated the observational data into three redshift bins 
and fitted, at each redshift, a linear relation between Rbr, 
Hhr, and the mass of the galaxies, in an identical fashion to 
that done by Azzollini et al. (2008). Using the linear fit, we 
measured the predicted values of Rhr and ^it>r for a galaxy 
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Figure 9. Colour profiles for several outputs of our simulation from z = 1 to z = 0. Dashed lines in each panel indicate the position of 
-Rbr m the V-band surface brightness profile (recall, Figure [4). Error bars reflect the RMS dispersion of the mean azimuthally-averaged 
colour in randomly-selected octants of the disc. The colour profiles show a minimum at the break radius, in agreement with observation. 



with the same mass as our disc at each of the different red- 
shift bins (note that the mass is different at each redshift) 
and their associated RMS dispersions. Figure UTI shows the 
predicted evolution of Rbr and fibr in the B-band for our sim- 
ulated disc, compared to the observed data. It can be seen 
that the evolution of these parameters within our simulated 
disc are very similar to that which is observed. 



4.6 Age and metallicity distribution in the 
outskirts 

Studies of stellar population in the outskirts of discs have 
found stellar populations which are of old-to-intermediate 
ages, with metallicity distributions peaking at relatively high 
metallicities ([Fe/H]~ — 0.7) (e.g., Ferguson & Johnson 
2001 (M31); Davidge 2003 (NGC 2403, M33); Galleti et al. 
2004 (M33)). It has been suggested that these properties 
are not expected in CDM models (e.g., Ferguson & Johnson 
2001), where the disc grows inside-out and at relatively re- 
cent epochs (z<l). Figure [T2"l shows the age and metallicity 
distributions of stars with Rbr < r < 15 kpc in our simulated 
disc. There is a very large scatter in the ages of the stars 
and a significant fraction of the stellar population has ages 
in excess of ~8 Gyr. The peak of the metallicity distribu- 
tion is near [Fe/H]~— 0.5 - i.e., even larger metallicities than 
that normally observed in the outskirts of galaxies. However, 
it should be noted that the "observed" metallicities are in- 
ferred via the use of colour-magnitude diagrams. When deal- 
ing with photometric metallicities, one must bear in mind 
that the underlying age distribution may affect the derived 
metallicity distribution. All of the studies cited above ob- 
tained their metallicites by comparing with the sequences 
of old Milky Way globular clusters. If a scatter in age is al- 
lowed, with a larger fraction of young stars (as we obtain in 
our simulation), the metallicity distribution will be skewed 
towards higher values and more closely resemble our derived 
distribution. We note that our metallicity distribution has a 
metal-poor tail, as reported in some empirical studies (e.g., 
Galleti et al. 2004). 



5 BREAK FORMATION 

As foreshadowed in Section 1, there is no clear consensus re- 
garding the underlying physical mechanism responsible for 
the surface brightness profile breaks observed both locally 
and at high-redshift. Broadly, there are two categories of 
theories - those related with the angular momentum dis- 
tribution of the stars and those related to star formation 
thresholds. 

In our simulated disc, the break is seen in the light 
profile, but not in the stellar mass density profile. This in- 
dicates (or at least suggests) that the break in our disc is 
likely related to a change in the stellar population prop- 
erties and, therefore, with a star formation threshold. The 
third and fourth rows of Figure U show the SFR density 
(with the SFR averaged over the last 1 Gyr) and gas sur- 
face density profiles, respectively. Examining Figure [4] we 
can see that the gas surface density at the break radius is 
not the same at all redshifts. For example, at z = 2.01 there 
is no clear break in the surface brightness profile, despite 
the gas surface density reaching values significantly below 
the star formation threshold. This would seem to contradict 
the aforementioned star for mation threshold sce nario as the 
driver behind the break. In lRoskar et al.l l|2008l ). the reason 
for the break was a sudden decrease in the gas density, but 
also, in agreement with our results, the density of the gas 
at that radius was higher than the star formation threshold 
imposed in their simulations. In their work, the decrease in 
the stellar density was due to an angular momentum limita- 
tion. By construction, in their simulations, the angular mo- 
mentum is directly proportional to the radius, which means 
that the high angular momentum material will take longer 
to cool. Therefore, the angular momentum determines the 
maximum extent of the gaseous disc and of the star-forming 
disc. In fully cosmological simulations, however, the accre- 
tion of matter is not as regular as the classical model of 
spherical shells collapsing from increasingly large radii. Par- 
ticularly, at high-redshift, before the galaxy is large enough 
to shock the accreting gas, accretion of gas occurs along fila- 
ments, as well as in mergers and clumps. The growth of our 
simulated disc is thus significantly more complicated than 
the disc in Roskar et al.'s idealised models. Although our disc 
does grow "inside-out", it does so in a less regular manner. 
In our simulations, in fact, the gaseous disc extends beyond 
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Figure 11. Predicted evolution with redshift of _Rb r and /^br (open circles) compared to the observed values from Azzolini et al. (2008) 
(filled circles with error bars). The error bars show the RMS-dispersion of the measured values for galaxies in a given redshift- and 
mass-bin. 



the stellar disc and, while in some timesteps the position of 
the break coincides with a change of slope in the surface gas 
den sity profile, this i s not true at all timesteps. 

iFerguson et al.l l|l998l ) proposed that the breaks in the 
star formation density could be due to an intrinsic correla- 
tion between azimuthally-averaged star formation rate and 
gas volume density, combined with a vertical flaring of the 
disc or a warp, in such a way that the transformation be- 
tween gas surface de nsity and volume d ensity varies with 
galactocentric radius (|Madore et al.|[T974l ). We plot, in Fig- 
ure [14] the ratio between the surface density and the volume 
density of the gas in our simulated galaxy. It can be seen 
that, indeed, there is a change in the slope of this ratio at 
(approximately) the break radius. 

Drops in the surface density of the gas at the onset of 
a warp have been obser ved in several galaxies (|J6zsal [20071 ; 
iGarcfa-Ruiz et all 12002), In fact, a possible connection be- 
tween truncated profiles and gaseous warps has been sug - 
gested several times in the literature (|van der Kruitll2007h . 
Figure [13] shows the edge-on projection of the gas density 
distribution of our disc at three different redshifts, match- 
ing those of Figure 3] The break radius coincides with the 
onset of a warp in the gaseous distribution in each of these 
three time-steps. Therefore, the decrease in the gaseous vol- 
ume density due to this flaring seems to be a fundamental 
mechanism responsible for the break seen in the average star 
formation rate and, hence, the break in the light profile. 

The coincidence of the initiation of the warp with the 
break radius at different timesteps is an encouraging piece 
of evidence in support of this suggestion. Having said that, 
this is but one simulation; a definitive statement awaits the 
analysis of a systemic suite of 20-30 comparable simulations, 
sampling a range of environments, assembly histories, and 
dynamical masses - such an ambitious program is currently 
underway within our group. 



6 RADIAL MIGRATION 

Figure [15] compares the final (at 2=0) galactocentric radii 
for stars in our simulated disc (Ranal) with the radius at 
which they were formed (Rinitiai). The right panel shows 
the same distribution, but only for those stars outside the 
break radius. Using disc stars with i?b r < r <15 kpc, we 



find that ~57% of the stars formed inside the break radius, 
~21% formed in situ, and ~22% formed at r >15 kpc or 
z > 3 kpc, including those formed elsewhere in external 
satellites - i.e., more than half of the stars currently in the 
outskirts of the disc formed at lower radii. 

The mean radial distance traversed by the disc star 
particles (all those with 3< r < 15 kpc and |z| < 3 kpc, 
discounting those formed in satellites for this calculation) 
is <|-Rfinai— J?initiai|>=l-7 kpc, while for those with R hr < 
r <15 kpc, this mean radial traversal distance is 3.4 kpc. 
It is clear that in our simulation there is considerable ra- 
dial migration of stars towards the external parts. The im- 
portance of such radial migration and mixing within the 
thin dis c has been the subject of sev e ral important stud- 
ies (e . g. ISpitzer fe Schwarzschild|[l953l; iBarbanis fc Woltiei 



.g. bpitzcr & be 
iFuchsl l200ll; IS 
20081 ; lHavwoodll2008l 



19671 ; iFuchs! l200ll; ISellwood fc Binnevl |2002| ; iRoskar et al 



It has been known for a number of years that scatter- 
ing by spiral structure and molecular clouds can gradually 
heat stellar discs, moving stars towards more inclined and 
eccentric orbits and changing the overall angular momen- 
tum distribution of the disc. Minor mergers and accretion of 
satellities can also produce heating (e.g., Velazquez & White 
1999, amongst many others). Stars that are on eccentric or- 
bits are at different radii at different phases of their orbits 
and, therefore, tend to naturally produce some radial mix- 
ing. However, the radial excursions due to this mechanism 
are not sufficient to explain the flat age-metallicity relation 
in the solar neighbourhood (Sellwood & Binney 2002). Typ- 
ical radial variations for a population of stars with a radial 
velocity dispersion ay are AR~ \^2a r /k, where k is the epy- 
ciclic frequency; for the old stars in the Milky Way the max- 
imum value for the excursions near the Sun is ~1-1.5 kpc. 

Another proposed mechanism for radial migration was 
called "churning" by Sellwood & Binney (2002), and is due 
to scatter ing of stars across co-rotation resonance by spi- 
ral waves ( Sell wood fc Binnevll2002l;ISellwood fc Pretolbool 
ISellwood fc Kosowskvll2002l ; lRoskar et al.ll200Sl ). resulting in 
a change in the orbit centres, but not the eccentricity. Such 
churning causes little increase in random motions (or heat- 
ing) because it preserves the overall distribution of angular 
momen tum (see Sellwood fc Binney 2002 for details). In the 
work of lRoskar et al.l l|2008l ). it was this churning mechanism 
that was predominantly responsible for the radial migration. 
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Figure 13. Edge-on projections of the logarithmic hydrogen density in H cm 3 , at three different redshifts (noted in the inset to each 
panel). Vertical arrows indicate the position of the break in the associated V-band stellar surface brightness profiles. 
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Figure 15. Final stellar galactocentric radius versus formation radius. The different panels correspond to different scalings on the 
abscissa. 



Our cosmological disc (as all hydrodynamical cosmolog- 
ical discs to date have been) is significantly hotter than that 
of the Milky Way (Gibson et al. 2008); as such, radial ex- 
cursion of stars in eccentric orbits is no doubt contributing 
to (and perhaps dominating) the mixing. In fact, the epyci- 
clic radius of the outer disc particles is 3.3 kpc, compatible 
with the mean calculated AR. We expect radial churning 
to have an important contribution too, although due to the 
larger scale-heigh of our disc, this c ontrib ution is certainly 
less important than in[Roskar et all i|2008T ). Furthemore, our 
cosmological disc is continuously being bombarded by small 
satellites. Younger et al. (2007) showed that merger-driven 
gas inflows deepen the central potential and contract the 
inner profile while, at the same time, angular momentum 
is transferred to large radius and causes the outer disc to 
expand. This has the net effect of driving radial migration 
of the stars towards the external parts. It is very difficult to 
disentangle all the different causes of radial migration in our 
cosmological disc(3 although it is certainly the case that each 
of these mechanisms is playing a role in both our simulated 
disc and in real galaxies. The main heating mechanisms in 
the disc presented here and in other cosmological simulations 
will be the subject of a future paper (House et al. 2009, in 
preparation) . 



Roskar ct al. (2008) showed that radial migration due 
to this churning mechanism was able to explain the ex- 
ponential decay of the surface brightness profile beyond 
i?br, the constancy of the break radius for stars of differ- 
ent ages , and predict age profiles in agreemen t with obser- 
vation s ( Bakos et aLlbOfiih . The simulation of iRoskar et al.l 
(2008) was performed under idealised conditions, with the 
disc growing by construction in an inside-out fashion, with 
very little star formation outside Rbi (which itself increases 
monotonically with time). Hence, the stars beyond the break 
in their model must have migrated there. While these con- 
trolled, high resolution simulations are necessary to study 
many of dynamical processes having place on the disc, it 
is also interesting (and important) to study the evolution 
of the surface brightness profile within a fully cosmological 
context, where hot and cold modes of gas accretion, and 
the effects of mergers and interactions are included. We find 
that the break radius does not always increase monotoni- 
cally with time (see, e.g., Figure ITT]) and the shape of the 
final profile is influenced by a variety of external processes. 
Therefore, in the next section, we explore the effect of radial 
migration on the overall properties of our simulated disc. 



5 Indeed, due to the large velocity dispersion of the stars com- 
pared with the Milky Way and to the larger scale-height, cosmo- 
logical simulations may not be the most appropriate for studying 
these particular mechanisms. 
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Figure 7. Surface brightness profiles measured in three different 
bands, as indicated in the inset to each panel. Dashed vertical 
lines show the position of the break radius, while the dashed- 
dotted lines represent the fits to the inner and outer exponentials. 



7 IMPLICATIONS OF RADIAL MIGRATION 
FOR DISC PROPERTIES 

7.1 Density profile 

What is the effect of the radial migration in the density pro- 
file of our simulated disc? Figure [16] shows the mass density 
profile at z — compared with the hypothetical density pro- 
file that we would observe if the stars had not migrated from 
their birth place. The most dramatic changes in the density 
profile happen in the outer parts. The movement of the stars 




R(kpc) 

Figure 8. Stellar surface density profiles for stars of different 
ages, as noted in the inset. The intensity of the break decreases 
with increasing age of the stellar population under consideration, 
with old stars showing an "upbend" at large radii. The total stel- 
lar surface density profile is consistent with a pure exponential 
out to ~5 exponential scale-lengths. 
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Figure 10. Mass-weighted age and metallicity gradients for our 
simulated disc at z=0. 



tends to weaken the intensity of the break that would have 
otherwise been seen in the stellar mass density profile. 

This behaviour is quite different to what is seen in the 
simulations where the redistribution of angular momentum 
is caused by a bar (eg., Foyle et al. 2008). Foyle et al. found 
that the outer profile in their simulations barely changed 
with time, while the inner profile became flatter due to the 
redist ribution of material. Nevertheless, IPohlen fc Truiillol 
(2006) distinguish two types of Type II profiles - those asso- 
ciated with the Outer Lindblad Resonance (OLR), and thus 
linked with the presence of a bar, and those not associated 
with bars, which are normally located at large radius (see 
Erwin et al. 2008 for a more detailed discussion) . It may well 
be that two different mechanism are responsible for these 
two variants of the Type II profiles. The profiles not associ- 
ated with the OLR are less common in early- ty pe galaxies 
than late-type galaxies |Pohren fc Truulloll2006l ) which may 
indicate that the presence of an extended gaseous disc are 
necessary to produce them. 

We showed, in Sec. 14.21 that the old stars have a density 
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Figure 12. Age and metallicity distributions for the stars with 
radius Rb r <^<15 kpc and |z| < 3 kpc. 
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Figure 14. Ratio between the surface mass density and the vol- 
ume density, as a function of galactocentric radius at redshift 
2=0, in our simulated disc galaxy. The position of the break in 
the surface brightness profile is indicated with a dashed line. 



profile with an upturn in the outer parts. Figure [T7] demon- 
strates that this is partially due to the significant effect that 
migration has on this outer profile. In this figure, the stellar 
density profile using different age bins is compared with the 
same profiles in the hypothetical situation where the stars 
have not migrated from the place of their birth. It is read- 
ily apparent that the surface density profile, in the absence 
of migration, appears truncated in all cases. This contrasts 
with the mass density profiles that we derived from the fi- 
nal output of our simulation, which appear much closer to a 
pure exponential (to at least ~5 exponential scale lengths). 
We conclude, therefore, that the exponential mass den- 
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Figure 16. Stellar surface brightness density profile of the disc 
at z = (dot-dashed curve) compared with the profile the disc 
would possess in the hypothetical situation that the redistribution 
of stellar material by secular processes was absent. 
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Figure 17. As in Figure 1161 but subdividing the stars by age, 
as indicated in the inset to each panel; upper panel: hypothetical 
situation of no radial migration; lower panel: as measured in our 
cosmological simulation. 



sity profile in Figure [4] would appear truncated/broken if 
there had not been radial migration. As old stars have had 
more time to migrate towards the external parts, the den- 
sity profile for the oldest age bin appears anti-truncated. 
This mechanism could also be responsible for the upbend- 
ing profiles seen in galaxies where a truncation in the star 
formation density does not occur (Type III profiles, recall- 
ing the nomenclature of Section 1). This was also proposed 
by Younger et al. (2007) who found that in their simula- 
tions, the radial migration was produced by angular momen- 
tum transfer during minor mergers; such a scenario receives 
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support from the empirical evidence associating asymme- 
tries and distortions with anti-truncated discs (Erwin 2005; 
Pohlen & Trujillo 2006). As mentioned in Sectiorj6] several 
mechanisms are certainly operating in our disc, in order to 
produce the radial migration of stars toward the external 
parts, including this effect of satellite accretion. 



7.2 Stellar populations 

Stars and interstellar gas in galaxies exhibit diverse chemical 
element abundance patterns that are shaped by their envi- 
ronment and formation histories. A wealth of surveys and 
satellite missions are devoted to obtain chemical patterns for 
individual stars as well as ages and kinematics in order to 
derive the star formation history and merger history of our 
galaxy (e.g., HIPPARCOS, RAVE, SEGUE, GAIA). The 
comparison with chemical and chemodynamical evolution 
models is expected to provide insights into the formation 
epoch of the different Galactic components and the relation 
between them. 

Chemical evolution models of the Milky Way usu- 
ally divide the "Galaxy" into ann uli with no radial trans- 
fer of material betwe en them (e.g. iFenner fc Gibson] 120031 ; 
IChiappini et al.l |2003T > . However, recent works (e.g. Roskar 
et al. 2008b; Schonrich & Binney 2008; Haywood 2008) 
has pointed out the importance that migration mechanisms 
might have in the studies of chemical evolution, in particu- 
lar, upon the observational constrains used to calibrated the 
models. For example, one of the most difficult observational 
results to reproduce with current "semi-numerical" chemi- 
cal evolution models is the lack of an apparent correlation 
between the age and the metallicity of the stars in the solar 
neighbourhood, as well as the large scatter in metallicity at 
a given age. On the other hand, full chemodynamical mod- 
els indicate that radial mi gration might be an explanation 
to this lack of correlation. IRoskar et alj (|200&f ) have shown 
that radial migration can largely affect the age and metallic- 
ity gradients, as well as the dearth of metal-poor stars in the 
solar neighbourhood (the so-called "G-dwarf problem"), ft 
has even been suggested that radial migration might be re- 
sponsible for the formation of the thick disc (e.g. Schonrich 
& Binney 2008, Haywood 2008). 

These studies provide the motivation to re-examine the 
stellar population distribution within our simulated disc, 
where infall of material and the accretion of satellites are 
taken into account naturally within our cosmological frame- 
work. 



7.2.1 The age-metallicity relation 

Chemical evolution models predict an increase of the metal 
content in the ISM as stellar generations die and pollute 
their surroundings with the byproducts of nucleosynthe- 
sis. One might expect therefore that the youngest (old- 
est) stars would also be the most metal-rich (metal-poor). 
However, there is little evidence of su ch an age-metallicity 
relation in the solar n eighbourhood IjFeltzing et alj l200ll ; 
iNordstrom et all |2004T ). Furthermore, the scatter in the 
relation is very large - i.e., at a given age, there is a 
large spread in the metallicity distribution. Classically, this 
was interpreted as a result of inefficient mixing of stellar 
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Figure 18. Age-metallicity relation for stellar particles with ra- 
dius between 7 and 9 kpc and scale-height \z\ < 3 kpc, chosen 
to represent a region roughly corresponding to the "solar neigh- 
bourhood". The red points indicate the values derived for the 
stars born in this region. 

ejecta. However, several authors (e.g., IWielen et al. I Il996l ; 
ISellwood fc Binnevll2002r ) have suggested the possibility that 
radial migration may have enabled old stars, formed at 
small galactocentric radii, to appear in the solar neighbor- 
hood, flattening any correlation between age and metallicity. 
IRoskar et ail (|2008r ) and Schonrich & Binney (2008) have 
shown that a flatter age-metallicity relation and a large 
spread in metallicities at a given age can be obtained with 
dynamical models where radial mixing is taken into account. 
Haywood (2006, 2008) observed that stars in the metal-poor 
and metal-rich end of the thin disc have orbital parameters 
which are offset from the main population, which he inter- 
preted as being the consequence of radial migration. 

Fig. 1181 shows the age-metallicity relation in the "solar 
neighborhood" (7 < r <9 kpc and z < 3 kpc) of our sim- 
ulated disc, compared with the age-metallicity relation that 
we would have measured in the absence of migration. As 
can be seen, our fully-cosmological study is in agreement 
with the aforementioned idealised models. Radial migration 
produces a considerable flattening, and an increase of the 
scatter, in the age-metallicity relation. In the absence of ra- 
dial migration, the stars in the solar neighborhood of our 
disc would show a relation between these two parameters. 

7.2.2 Age gradient 

We showed, in Sec. 14.41 that the radial age profile has a 
minimum at the b reak radius, similar to that found by 
IRoskar et alj ((2008) . In their work, the radial extent of the 
star forming disc was limited by the maximum angular mo- 
mentum of material that was able to cool at each timestep. 
Therefore, by construction, the disc grows inside-out and 
the youngest stars are situated at the break radius. The 
"up-bend" in the age profile at r>R br was produced, exclu- 
sively, by stars that had migrated from the internal parts of 
the disc. Because older stars have more time to travel larger 
distances, the trend of the age with radius reversed after the 
break radius. 

To explore if this is also true in our simulation, we plot, 
in Fig. 1191 the mass-weighted age profile that the stellar 
disc would have if the stars did not migrate from their birth 
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place, compared with that derived from the final at z = 0. 
As can be seen, even in the absence of radial migration, a 
'U-shape' profile in the age gradient is still visible. Stellar mi- 
gration changes the shape of this profile, but the increasing 
mass-weighted age with increasing radius in the outer parts 
is not due to radial migration. Figure shows the evolution 
with time of the total star formation rate density for stars 
formed at different galactocentric radii, including regions in- 
side and beyond the break radius. For radii 3<r<5 kpc, the 
SFR density decreases with time, while the opposite hap- 
pens for 5<r<7 kpc. The consequence of this is that the 
the fraction of young-to-old stars increases with radius until 
r=Rbr, consistent with expectations, as the early accreted 
mass has low specific angular momentum. These trends are 
responsible for the decreasing mean age with radius and are 
in agreement with the inside-out scenario for the formation 
of disc galaxies (Ryder & Dopita f994). However, the trend 
dissapears beyond the break radius, where the star forma- 
tion rate is low and essentially constant throughout the evo- 
lution of the galaxy. The difference in the star formation rate 
between the outer (f0< r < 15 kpc) and inner (7<r<9 kpc) 
parts is higher at later times, and that is the main reason for 
the upturn in the age gradient. Another way to appreciate 
this is by studying the radial profile of the birth parameter, 
defined as the cur rent versus averaged s tar formation rate 
fr=fSFR/<SFR>) (|Kennicutt et al]|l994 ). Instead of plot- 
ting the current star formation rate, we averaged the star 
formation over the last 1 Gyr. This is mainly due to the 
relative few number of particles with ages below this value 
outside the break radius. We plot, in Fig. 1211 this parame- 
ter as a function of galactocentric radius. As can be seen, 
b increases almost linearly with radius, as expected for an 
inside-out formation scenario, until the break radius, where 
it reaches a plateau and then decreases. The increase in the 
error bars at the break radius reflects the asymmetries in 
the age distribution of the stars beyond this radius (recall 
Fig. El. 
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Figure 19. Mass-weighted azimuthally-averaged stellar age and 
metallicity gradients. Solid lines: Theoretical gradient for the hy- 
pothetical case were stars do not migrate from their birth place. 
Dashed-dotted line: Profile measured in the final timestep of our 
simulation. Dotted line: Profile measured in the final timestep of 
the simulation after eliminating those stars which formed out- 
side the disc (those with initial galactocentric radii in excess of 
25 kpc). 



We argue that the U-shape age profile is the direct con- 
sequence of the existence of a break in the star formation 
density. If the star formation outside the break had not de- 
creased suddenly, the age gradient would decrease, or re- 
main constant, until th e edge of the optical disc. This is sup- 
ported by the result of IBakos et al.l (120081 ) who only found 
the tell-tale U-shaped colour profiles for galaxies possess- 
ing a Type II profile. The galaxies with an essentially pure 
exponential profile within their sample showed a plateau 
(and not an up-bend) in the colours at large radii. This 
up-bending age profile does not mean that the disc did not 
form inside-out. In fact, the "overall" formation of the disc 
remains inside-out. However, in our disc, the decrease of star 
formation in the external parts - due to a decrease in the 
volume density of the gas - results in redder colors beyond 
the break radius. 



10 <r < 15 kpc 




4 6 8 10 12 

t(Gyr) 

Figure 20. Comparison of the star formation rate density with 
time in 4 different regions of the galaxy disc - black: between 10 
and 15 kpc; dark red: between 7 and 9 kpc; purple: between 5 
and 7 kpc and orange between 3 and 5 kpc. 



In Fig. [I9]we also compare the age profile of the galaxy 
with the one it would have if all the stars formed in the disc 
- i.e., if satellites were not accreted. It is apparent that for 
this case, the accretion of satellites has little effect on the 
stellar population gradients. 
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R(kpc) 

Figure 21. Star formation over the last 1 Gyr divided by the 
total mass of stars formed before this epoch (the so-called "birth 
rate") as a function of galactocentric radius. 



7.2.3 Metallicity gradient 

The metallicity gradient in the disc and its evolution with 
time provides constraints to our understanding of the forma- 
tion and evolution of galaxies. The presence of a metallicity 
gradient in the Milky Way is widely accep ted, although its 
exact slope and shape remain contentious l|Chiappini et al.l 
l200ll : lAndrievskv et alJl2002h . A related issue is the evolu- 
tion of this gradient with time; this has been approached 
from both the theoretical (Gotz & Koppen 1992; Koppen 
1994; Molla et al. 1997; Henry & Worthey 1999; Chi- 
appini et al. 2001) and observational perspectives (Friel 
et al. 2002; Maciel 2001; Maciel, Costa & Uchida 2003; 
Stanghellini et al. 2006) but it is still not clear whether 
the metallicity gradient in our Galaxy flattens or steepens 
with time. Measurements using HII regions, B-stars, and 



planetary nebulae, find gradients ranging 


from ~— 0.04 to 


~— 0.07 dex kpc 1 jDeharven 


e et al.l 2000 


; Afficrbach et al. 


1997: Dafion & Cunhal 12004 


: ICummcrsbach et al.l !l998: 


Costa et al.l 20041: Maciel & Quirezal 19991. although see 



Stang hellini et al. 2006). Even more recently. ICarraro et al.l 
(2007), using five old open clusters, suggest an even shal- 
lower gradient of —0.018 dex kpc -1 . There is also contro- 
versy in the literature about the shape of the gradient, in- 
cluding the presence or not of a discontinuity at around 
10 — 12 kpc. Part of the uncertainty in this field can no 
doubt be traced to the necessary use of disparate tracers of 
the gradients, each which probe a different age (or range of 
ages). 

In our simulations, we find that the metallicity gradient 
is flatter for old stars than it is for young stars, as shown in 
Figure 1221 It can also be seen that radial migration is par- 
tially respo nsible for this flatte ning (see Fig. I19J1 , in agree- 
ment with iRoskar et alj (|2008l ). We also predict that the 
optimal galactocentric radius in which to seek differences 
between the gradients in the old versus young populations 
is 3<r<5 kpc. At can be seen in Figure [T9l (lower panel), the 
mean metallicity gradient in our simulation shows a flatten- 
ing at r>12 kpc. 
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Figure 22. Azimuthally-averaged metallicity profile at redshift 
z=0 for stars in different age bins, as noted in the inset. 



8 DISCUSSION AND CONCLUSIONS 

We have studied the origin of the disc truncations and the 
stellar population properties of a disc galaxy formed within 
a cosmological framework. To do this, we have analysed a 
mult i- resolved cosmological simulation of a Milky Way-mass 
halo, including dark matter and gas dynamics, metallicity- 
dependent cooling, UV heating, star formation, and super- 
novae feedback, using the AMR code RAMSES. Our simulated 
disc shows a break in the surface brightness profile at ~3 ex- 
ponential scale- lengths, similar to those observed in a large 
fraction of disc galaxies. The position of the break does not 
coincide with the radius at which the gas density reaches 
the canonical threshold for star formation. In our simula- 
tion, this break in the stellar light is due to a decrease in 
the star formation density per unit area (averaged over the 
past ~1 Gyr). This decrease in the star formation density 
originates from a decrease in the volume density of gas at the 
break radius, which itself coincides with the radius at which 
the gas disc begins to warp (R W ar P ). A relationship between 
truncations and warps has been pointed out by several stud- 
ies (see van der Kruit 2008 for a review as well as van der 
Kruit 2001). van der Kruit (2007) found that the distribu- 
tion of Rwarp/Rbr in a sample of SDSS galaxies was statis- 
tically consistent with all warps starting at ~1.1 Rbr- Our 
simulation is entirely consistent with these observations. At 
a first sight, therefore, it seems that the presence of a warp 
is a condition for the presence of a break in the stellar light 
distribution. However, we re-iterate that this need to be con- 
firmed with a larger sample of simulated and observed discs. 

We analysed the redistribution of material and angular 
momentum in the disc, finding that it affects considerably 
the final shape of the mass density profile, especially in the 
outskirts. In fact, 57% of the stars with r> i?br formed at 
lower radii. The reason why the surface density profile does 
not show a truncation similar to the one observed in the stel- 
lar light is traced to this migration of stars towards the outer 
disc. We suggest that truncated and anti-truncated profiles 
can be produced by a different combination of two processes: 
(1) a change in the slope of star formation profile with ra- 
dius due to a change of the slope in the gas density profile 
(which causes the truncation in the light) , and (2) radial mi- 
gration of stars formed in the internal parts towards radii in 
excess of Rbr- We speculate that up-bending profiles can be 
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produced in our simulation when the distribution of the gas 
changes smoothly with radius. In that case, process (1) is 
suppressed, the star formation rate per unit area in the disc 
changes smoothly with radius, and the only deviations from 
a pure exponential profile are due to stars migrating from 
disc interior, producing an increase in both the surface light 
and mass density in the external parts. This would make the 
external parts somewhat redder, in agreement with observa- 
tions (Bakos et al. 2008). Pure exponential profiles can be 
similar to anti-truncated profiles, but with less migration. 

We have also studied the stellar populations in the disc 
and the influence of the migration of stars and the accretion 
of satellites in the age and metallicity profiles. We found, in 
agreement with recent observations (Bakos et al. 2008) that 
the disc shows a U-shaped age profile, reaching the minimum 
age at the position of the break. It has been proposed in re- 
cent works that this profile is produced by the combination 
of two process: (1) the inside-out growth of the disc (result- 
ing in the age decreasing with radius until Rbr), and (2) the 
migration of old stars from the interior to the exterior parts 
of the disc beyond Rbr- In this picture, the increase of age 
with radius beyond the break results from older stars hav- 
ing more time to travel from the inner disc and, therefore, 
reach a greater galactocentric radius. These models predict 
a flat metallicity gradient for the old stellar populations, as 
these stars are mixed very efficiently. In our simulation, the 
U-shaped profile is due to a different rate of star formation 
between the internal and the external regions (inside and 
outside the break radius) . Migration of the stars also has an 
impact on this profile, but the U-shape age profile with the 
minimum at the break radius appears even when migration 
is suppressed in the disc. We find that the U-shaped age pro- 
file is due to the same mechanism as the one producing the 
break in the surface brightness distribution - ie, a diminu- 
tion of the star formation rate in the external region, due 
to the warping of the gaseous disc and, therefore we predict 
that this type of color profile is only present when there is 
a break in the light distribution. This is in agreement with 
the recent work of Bakos et al. (2008), at low redshifts, and 
Azzolini et al. (2008) at redshifts out to z ~1.1. 

The distribution of ages in the outer disc of our simu- 
lation is very broad with stars nearly as old as the age of 
the Universe itself. The presence of old stars (>I0 Gyr) in 
the outer disc of some nearby galaxies, such as M31 (Fergu- 
son & Jonhson 2001), has been used as an argument against 
ACDM models, arguing that, if feedback is included in order 
to produce large discs, the formation of the disc is delayed 
and the resulting stellar populations should be of young and 
intermediate ages. We show here that in a ACDM frame- 
work we were able to produce disc with both a realistic 
scale-length and a considerable number of old stars in its 
external parts. The metallicity distribution of stars in the 
outskirts of the disc peaks at [Fe/H]~ —0.5 and shows also 
a significant spread, also consistent with observations. 

Finally, we studied the age-metallicity relationship in 
a representative "solar neighbourhood" of our simulation. 
There have been claims suggesting that radial mixing of 
stars might be responsible for the absence of an obvious 
correlation between these two parameters in the solar neigh- 
bourhood, and for the large scatter in metallicity at a given 
age. We agree with these claims, showing that a flat rela- 
tionship between age and metallicity, and significant scatter 



in metallicity at a given age, are a natural outcome within 
our cosmological simulation. 

In a forthcoming paper we will study the origin of the 
warp in the gaseous distribution, its lifetime, asymmetry, 
and intensity. This will be supplemented with a suite of 20- 
30 additional disc simulations, spanning a range of mass, 
environment, and assembly history. This will allow us to 
provide significantly more robust predictions for the rela- 
tionship between break radii and other empirical character- 
istics of galaxies. 
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APPENDIX A: DEPENDENCE OF THE BREAK 
RADIUS UPON INCLINATION 

Recently, van der Kruit (2008) studied the correlations be- 
tween break radius and other galaxy properties, finding fun- 
damental differences between the trends described by face- 
on and edge-on galaxies. In particular, he did not find a cor- 
relation between the position of the break and the rotational 
velocity in the face-on sample, while the correlation was ev- 
ident in the edge-on galaxies. In principle, breaks are more 
easily detected in edge-on galaxies, due to their favored in- 
tegrate line-of-sight vantage point. Face-on breaks are more 
difficult to detect and if the truncation radius changes with 
azimuth, the azimuthally-averaged profiles would "smear 
out" the intensity of the break. Furthermore, intrinsic de- 
viations from circular symmetry, such as spiral arms (espe- 
cially for younger populations) should complicate matters 
for less-inclined systems. On the other hand, edge-on stud- 
ies are more susceptible to line-of-sight effects caused by 
integrating through the disc. The differences between the 
parameters measured in face-on and edge-on samples can 
be obtained statistically, but obviously, they are difficult to 
quantify. In this appendix, we compare the break radius and 
the intensity of the break in different bands by viewing our 
simulation various face-on and edge-on configurations, in or- 
der to quantify possible differences due to geometry. Differ- 
ences due to the varying dust contribution with inclination 
cannot be treated here, but we refer the reader to Pohlen et 
al. (2008) for a thorough analysis of this issue. 

All the profiles previously showed in this work have 
been obtained by azimuthally-averaging the luminosities of 
the stars within 3 kpc of the disc mid-plane. Observers, of 
course, integrate their luminosities along the entire line-of- 
sight. In order to mimic this we compare the profiles inte- 
grating the light of all stars within a distance of 10 kpc from 
the disc mid-plane. Figur jAll shows the surface brightness 
profiles in different bands obtained for face-on and edge-on 
orientations. Table IA1I shows the measured values for the 
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Figure Al. Surface brightness profiles in different bands cal- 
culated in face-on (azimuthally averaging) and edge-on configu- 
rations. For the face-on projection we have integrated between 
different distances from the disc mid-plane, as indicated in the 
inset. 

position of the break and the angle between the two expo- 
nential fits in each configuration. 

We do not find any fundamental difference between the 
position of the break or its intensity for the face-on and 
edge-on values; they agree, within the uncertainties. The ta- 
ble also shows the measured inner and outer scale- lengths. 
The uncertainties in the measurement of the scale-length 
are <0.1% and are, therefore, not listed. We can see that 
contrary to what happens with the position of the break, 
we did find differences in the scale-lengths of the disc when 
measured at different inclinations. In particular, the scale- 
lengths are systematically lower in the edge-on projections, 
counter to that claimed by van der Kruit (2008). Due to the 
smaller scale-length obtained in edge-on systems, the RsR/h 
measured for this projection is larger than for the face-on 
projection. These differences have been found previously in 
empirical studies, but it is not clear if they are due to pro- 
jection effects or to different techniques to mark the break 
(see Trujillo & Pohlen 2006). 

Finally, we do not find a sharper edge in edge-on galax- 
ies than for face-on systems, for which we conclude that this 
effect, reported in observations, must be the consequence of 
dust extinction. 
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band 


face-on 


face-on 


edge-on 






|z <10 kpc 


|z < 3 kpc 




Rbr 


B 


9.3 ±0.3 


9.6 ±0.7 


9.8 ± 1.0 


^inn 


13 


2.83 


2.75 


1.95 




B 


1.74 


1.46 


1.30 


angle 


B 


10.2 


14.4 


10.4 


R br 


V 


9.2 ±0.4 


9.5 ±0.6 


9.3 ± 1.3 




V 


2.61 


2.53 


1.81 




V 


1.83 


1.47 


1.41 


angle 


V 


7.6 


12.5 


6.3 


Rbr 


R 


8.9 ±0.7 


9.8 ±0.4 


9.2 ± 1.7 


angle 


R 


5.9 


11.1 


6.5 


hznn 
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2.51 


2.43 


1.83 


hout 


R 


1.90 


1.52 


1.42 


Rbr 
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8.7 ± 1.1 


9.6 ±0.4 


8.3 ±1.0 


angle 
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4.6 


9.3 


4.2 
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2.42 


2.34 


1.80 
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1.96 


1.57 


1.52 


Rbr 
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8.2 ±1.7 


9.5 ±0.3 


7.8 ± 1.4 


angle 
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3.4 


8.2 


3.4 
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2.31 


2.24 


1.76 


hout 
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1.99 


1.58 


1.54 


Rbr 
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8.1 ±2.9 


9.5 ±0.3 


12.1 ± 6.0 
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2.20 


2.20 


1.60 
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1.59 


1.52 


angle 
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9.2 ±2.1 


9.4 ±0.3 
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2.19 


1.65 


hout 
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1.95 


1.58 


1.55 


angle 
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2.8 


7.7 






Table Al. Break radius and angle between the two exponential 
fits in different bands. 



